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Abstract. We introduce a new method to measure frequency separations and mode hfetimes of stochastically excited 
and damped oscillati ons, so-called solar- like oscillations. Our method shows that velocity data of the red giant star ^Hya 
iprandsen et alJl2002h support a large frequency separation between modes of roughly 7 //Hz. We also conclude that the data 
are consistent with a mode lifetime of 2 days, which is so short relative to its pulsation period that none of the observed fre- 
quencies are unambiguous. Hence, we argue that the maximum asteroseismic output that can be obtained from these data is an 
average large frequency separation, the oscillation amplitude and the average mode lifetime. However, the significant discrep- 
ancy between the theoretical calculations of the mode lifetime iHoudek & Gou gh 2002) and our result based on the observations 
of ^Hya, implies that red giant stars can help us better understand the damping and driving mechanisms of solar-like p-modes 
by convection. 
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1. Introduction 

The mode lifetime of solar-like oscillations is an important pa- 
rameter The interpretation of the measured oscillation frequen- 
■ cies (and their scatter) relies very much on knowing the mode 
lifetime, but currently we know very little about how this prop- 
erty depends on the stellar parameters (mass, age and chemical 
composition). The theoretical estimates of mode lifetimes are 
based on a simplified description of the convective environment 
in which the damping and excitation of the modes takes place. 
Measurements of the mode lifetime, t, in different stars will 
be very helpful for a more thorough treatment of convection in 
stellar modeling. In this paper the mode lifetime refers to the 
time for the amplitude to decrease by a factor of e. 

The number of measurements of the mode lifetime, or 
damping time, is still very limited. Observations of main- 
sequence stars imply mode lifetimes of a few d ays in the Sun 
/Libbrecht'1988; 'Chaplin et al."l99f), ceCen A jBedding et all 
[2004) and a Cen B (Kieldsenet al. 2005). As pointed out by 
bedding et al] ll2005hl) . independent observational studies on 
the star Procyon do not agree on the measured frequencies, 
a disagreement that could be the result of a short mode life- 
time. The power sp ectrum of the K giant Arcturus reported by 
iRetter et'al] ( l2003h could also be explained as a short mode 
lifetime (t - 2 days) of a single mode. If the mode lifetime 
does not increase with oscillation period, this would limit the 
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prospects of asteroseismology on evolved stars that have peri- 
ods of several hours or longer, because of poor coherence of 
the oscillations. Only when we look at M giants - the semi- 
regular variables - do we see evidence for longer lifeti mes, 
ranging from years to decades fchristensen-Dalsgaard e t alJ 
l200lUBedding et alJ2003ir2005ah . 

The theoretical predictions of mode lifetimes for un- 
evolved stars like the Sun (Houdeket al. 1999) and a Cen A 
dSamadi et alJ l2004) are in a fairly good agreement with the 
observed values. However, for the more evolved red giant star 
^Hya, there seems to be a significant discrepancy between the- 
ory (r ~ 15-20 davs:iHou dek & Gough 2002) and observation 
(r ~ 2-3 davs: Istello et al.i .2004. hereafter Paper I). In Paper 

1 we also measured the oscillation amplitude to be roughly 

2 ms~', which was in goo d agreement with the theoretical value 
llHoudek & Goughl2002h . 

In Fig. [2 we plot the measured ratios between the mode 
lifetime and period (the oscillation 'quality' factor) as a func- 
tion of period for selected stars, including the theoretical value 
for ^ Hya. Roughly speaking, the 'quality' factor is the number 
of oscillation cycles with constant phase, and the higher this 
number, the better we can determine the frequency. Note that 
the relation between the mode lifetime r and the the FWHM 
r (in cyclic frequency) of the corresponding resonant peak is 
r= l/(7rT). 

In this paper we further investigate the mode lifetime of 
^ Hya and examine whether it limits the possible astrophysical 
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Fig. 1. The oscillation 'quality' factor vs. period for se- 
lected stars. Filled symbols are measured, while the empty 
symbol shows the theoretical value. The arrow indicates a 
lower limit. Luminosity classes are indicated at th e top. 
The measured values are from: o-Cen B, Kieldsen^^^ 



1 20051) : Sun, IChapUn et all (Il997l): aCen A, Kieldsen et al 



1 2005h : ^Hya ftheoreticaP. iHoudek & GoueHI (.200Z) : ^Hya 
(observed). Paper I; Arcturus. iRetter et" [J ll2003h: L 7 Pup, 
iBeddinget aDoOOSal) : SV Lyn and R Dor. lDindl ( l2004 . 



output. We introduce a new method that measures repeated fre- 
quency patterns (e.g. the large frequency separation, Avq) and 
the mode lifetime of solar-like p-mode oscillations. We further- 
more analyze the stellar power spectrum to establish the most 
likely mode lifetime, by comparing the cumulative power dis- 
tribution of the observations with simulations. Finally, the am- 
biguity of the measured frequencies is quantified. 

2. The ^ Hya data set and previous results 

We use the same data set as in Paper I, which comprises 
433 measurements of radial velocity covering almost 30 days 
of single-site observations using the Coralie spectrograph at 
La Silla (ESO, Chile). The average noise per measurement is 
o"measure = 2.33 ms"', and the data provide a clear detection 
of excess power in the Fourier spectrum at roughly 90juHz 
(period ~ 3 hours) with many peaks of S/N> 3 (see Fig. Et. 
For fu rther details about the observations see iFrandsen et alJ 
ll2002h . Note that during t he red uction the data was high- 
pass filtered (Frandsen et al.] Eo02l) . but we have shown with 
simulations that this does not affect the power spectrum 
abov e 20AiHz. Us ing the autocorrelation of the power spec- 
trum iFrandsen et al. (2002) found a frequency spacing of 
6.8 //Hz, which is in good agreement with the large fre- 
quency separation, Avo, from a pulsation model of the star 
(Christensen-Dalsgaard 200-^. F or further ana l ysis on the stel- 
lar p arameters of f Hya see e.g. Stello (2002); Frandsen et al. 
( E00Z) : .Teixeii-a et aL (.2003 ): Paper h Thevenin et al.. (.2005) . 

In Table Q] we give the frequencies used in the current in- 
vestigation. We used the conventional method of iterative sine- 
wave fitting Cprewhitening') to measure the 10 frequencies 
listed in Table ITI (se e Appendix | A|. The uncertain t ies are cal- 
culated according to Montgomerv & 0'Donoghuel(ll999l) (see 
also .Kieldsea.2003) . We note that these frequencies are n ot ex- 
actly the same as those quoted bv IFrandsen et alJ (l2002l) . who 
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Fig. 2. Power spectrum of ^Hya. The frequency axis of the 
spectral window in the inset is scaled to match that of the main 
plot. 

Table 1. Measured frequencies and ampUtudes for ^ Hya 



ID 


Frequency 
pHz 


Amplitude 

ms"' 


S/N 


Frandsen et al. 
2002 


Vl 


85.91(3)" 


1.89(24)" 


6.3 


85.96(3) 


V2 


59.43(3) 


1.75 


5.8 


59.43(3) 


V3 


79.13(3) 


1.65 


5.5 


79.13(3) 


V4 


95.21(4) 


1.33 


4.4 


95.28(4) 


V5 


108.20(4) 


1.24 


4.1 


108.22(4) 


V6 


101.32(4) 


1.17 


3.9 


101.16(4) 


V7 


98.76(4) 


1.12 


3.7 


98.77(5) 


V8 


112.29(5) 


1.11 


3.7 




V, 


105.14(5) 


1.01 


3.4 


105.13(5) 


VlO 


73.38(5) 


1.02 


3.4 




^ The errors are derived from: o"/ = ^/6/(T^,^y^■S/N■n^^^), 
and o-a„,p = V277r</imp), where </J,mp> = o"mcasure Vir/iVobs, 



and N„t,s = 433 (see text). 



used a different method to measure the frequencies. This dis- 
crepancy is not surprising because the alias peaks (Fig.|3 are 
Ukely to introduce small frequency differences, depending on 
the method used to extract the frequencies. The uncertainties 
indicated in Tableware based on the signal-to-noise (S/N) and 
do not take into account the scatter of the frequencies arising 
from a short mode lifetime. This extra frequency scatter is the 
main subject of this paper 



3. The method 

3.1. Background and definitions 

From theoretical calculations we expect the oscillations 
in red giant stars to be do mi nated by radia l modes 
(IChristens en-DaIsg , aardl l2004 lOuenther et all l200(A 
IPziembow skietalJ l200lh . which is also supported by 
the o bservations o f f Hya (Paper I) and other red giant s 
jBuzasi et aL.200a: iKiss & Beddindl2003t iRetter et alJl2003l) . 
We therefore assume that the frequencies, v„, of ^Hya wiU 
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Fig. 3. Schematic illustration of the frequency scatter due to a 
finite mode lifetime. Dashed lines are the true mode frequen- 
cies (Eq.[0, and solid peaks are measured frequencies. 

show a simple comb pattern in the power spectrum with only 
radial modes. This can be described by the linear relation 



AvqW + Xq , 



(1) 



where Avq is the separation between mode frequencies of suc- 
cessive order n and Xq is an offset. 

The finite lifetimes of the oscillation modes will intro- 
duce deviations of the measured frequencies from the true 
mode frequencies and hence from the reg ular comb pattern 
jAnderson et alJ ll990': 'Bed ding et al.ir2004 . In Fig. I^we il- 
lustrate this by showing, schematically, the comb pattern of 
'true' mode frequencies (dashed lines; Eq.Q and the measured 
frequencies (solid peaks). The deviations, indicated with X(v,) 
(Eq. |2}, are independent because each mode is excited inde- 
pendently. Shorterjnode lifetim es \ yill give larger deviati ons. 
Following Bedding et al. (2004) and Kieldsen et al. ( 2005), we 
estimate the mode lifetime from the scatter of the measured fre- 
quencies about a regular pattern. Our method differs by using 
a comb pattern for the reference frequencies, and by not re- 
quiring us to assign the mode order or degree to the measured 
frequencies. This is an advantage when the power spectrum is 
crowded due to aliasing and noise, as in this case, but the draw- 
back is lower sensitivity to the mode lifetime. 

The first step is to find the comb pattern that best matches 
our measured frequencies, v = [vi , . . . , vn]. We do this by min- 
imizing the RMS difference between the measured frequencies 
and those given by Eq.^ with Avq and Xq as free parameters. 
The individual deviations are 



X(vi) = (v; - Xo) mod Av() , 



(2) 



where we use a modulo operator that returns values between 
-jAvo and iAvo. The RMS scatter of X is then 



V2 

crx(v;Avo,Xo) = - — 
Avo \ 



/=1 



(3) 



Here, = 10 is the number of frequencies. The normalization 
factor, V2/AV0, is included to make ctx independent of Avq for 
randomly distributed frequencies. We find the mini mum of cry, 
min(crx), using the Amoeba minimization method dPress et alJ 



Fig. 4. 
(Eq.E). 



Deviation of input frequencies from a comb pattern 



Il992h . and calibrate it against simulations with known mode 
lifetime, as described in the next section. As a by-product, we 
also get estimates for Avo and Xq. 

3.2. Simulations and calibration 

We simulated the ^Hya time series using the method described 
in Paper I. The oscillation mode lifetime was an adjustable pa- 
rameter, assumed to be independent of frequency, while the 
other inputs for the simulator were fixed and chosen to repro- 
duce the observations (see Paper I, Fig. 12). To make the sim- 
ulations as realistic as possible, the input freque ncies were the 
radial modes from a pulsatio n model of ^Hya (iiTeixeira et alJ 
I2003tlchristensen-Dalsgaardi 2004). Since our model assumes 
the frequencies to be strictly regular (Eq.O, the intrinsic devi- 
ation of the input frequencies from a comb pattern (see Fig.|3 
will contribute to ax- However, as described in Sect. 13.31 this 
contribution turns out to be negligible and therefore does not 
affect our measurement of the mode lifetime in ^ Hya. 

For different values of the mode lifetime, we simulated 100 
time series with different random number seeds. For each we 
measured 10 frequencies (vi, . . . , vio) using iterative sine-wave 
fitting (see Appendix lAb and then minimized crx- This pro- 
vided 100 values of min(crx) for each mode lifetime, which 
can be compared with the observations. In Figs.|5]and|6lwe il- 
lustrate the method, and Figs. and |8l show the results. Fig.|5] 
shows X{vi) of all 1000 frequencies (10 frequencies from each 
of 100 simulated time series) for two mode Ufetimes, (a) r = 17 
days, and (b) t = 2 days. The difference in mode lifetime is 
clearly reflected in the difference in the scatter of X(v,). We 
note that the distribution of X(v/) (right panels) is affected by 
the presence of false detections of alias peaks from neighbour- 
ing modes at roughly X(vi) ~ 2 juHz and -2 fjHz (see Fig.|5t). 

In Fig. |6l we show crxiv; Avo,X^)) for the case of t = 
17 days. To smooth the plotted surface we show the aver- 
age of all 100 simulations. However, during the minimiza- 
tion of (Tx, each simulations was treated separately. Note that 
crx(v; Avo,Xq) = crxiv; Avo,X() + nAvo), where n is an integer, 
and that crx has maxima and minima for roughly the same Avq 
that are separated by ^Avq on the Xo-axis. 
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Fig. 5. The deviation from a comb pattern (Eq.|2} for 10 mea- 
sured frequencies from each of 100 simulated time series. The 
two horizontal bands (at X(v,) ~ 2 fiHz and -2 fiHz, most evi- 
dent in panel (a)) below and above the central band are signa- 
tures of the alias peaks that are present in our spectral window 
(see Fig.|2jl. To the right, the distributions are shown for each 
panel. 



The output parameters min(crx) and Avo are plotted in 
Fig.0 We see clearly that a smaller min(crjf) (more prominent 
comb pattern) gives a more accurate Avq determination. In the 
case of short mode lifetime, the frequency pattern is generally 
not very pronounced (high mm(crx)) mainly due to the false 
detections from alias peaks. 

Due to the stochastic variations in the simulated time series, 
the value of min(crx) (Fig. shows a large intrinsic scatter 
from one simulation to the next. This tells us the precision with 
which we can determine the mode lifetime, r, from a single 
data set. For each r, we compared the observed min(cr;f) with 
the distribution from 100 independent simulations. This gives 
a measure of whether the observed value is consistent with t. 

Fig. |8l shows the distributions of mm(crx) for simula- 
tions as histogram-plots, for different mode lifetimes (r = 
17 days, 5 days, 2 days, and 6 hours). The dashed line indicates 
the observed value. For a mode lifetime of 17 days (corre- 
sponding to the damping rate in cyclic frequency, plotted by 
lHoudek & Gough 200^ of 0.1 fMz) only a few out of 100 tri- 
als have a value for min(cr;i^) as high as the observations. 

The exact number depends on the exact choice of analysis 
method (see Sect. 13. 3t . We find the best match with observa- 
tions for mode lifetimes of about 2 days, in good agreement 
with Paper I. However, we note that there is also a reasonable 
match for all mode lifetimes less than a day, as their distribu- 
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Fig. 6. Top panel: 3D surface of crx, (averaged from 100 sim- 
ulations with T = 17 days) as a function of the free parame- 
ters Xo (plot range 0,20) and Avo (plot range 5.5,8.0). Middle 
panel: The surface of crx viewed from above, where contours 
are plotted to guide the eye. Solid lines indicate contours above 
the mean level, while the dotted lines are those below. Bottom 
panel: The outline of the 3D surface viewed edge on in the 
Xo-direction. The two output parameters of the minimization 
process, min(crx) and the corresponding Avq, are indicated by 
dashed lines. 



tions all look very similar to the bottom panel. Randomly dis- 
tributed peaks also show similar distributions to t = 6 hours. 
Hence, if the mode lifetime of ^ Hya is only a fraction of a day 
it would definitely destroy any prospects for asteroseismology 
on this star. We have confirmed this with simulations that do not 
include any comb patterns by simulating a single mode with a 
mode lifetime of 2 hours, similar noise, excess power hump, 
and with the same total power as the above simulations. 
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Fig. 7. Output from our method applied on 100 simulations of 
^Hya with mode lifetime 2 days (diamonds), and 17 days (plus 
symbols). The square shows the observational point, and the 
size of the symbol indicates the uncertainty according to the 
chosen CLEAN algorithm (cf. Sect.|33J. 

3.3. Robustness of method 

An important step for a correct interpretation of the results 
shown in Fig. |8]is to test the robustness of our method. We 
tested the dependence of the two relevant output parameters 
Av() and min(crx) on the following: 

1 . the initial guesses of Avq and Xq in the minimization pro- 
cess of CTx, 

2. the number of frequencies used to calculate crx, 

3. whether weighting of frequencies is used when calculating 
crx, 

4. the frequency separation of input frequencies, 

5. the intrinsic scatter of input frequencies, 

6. the method for measuring frequencies. 

We discuss each of these in turn in AppendixlEI but in sum- 
mary, none of the listed points (1-6) have a significant effect on 
the results shown in Fig.|8l 

4. Power and frequency analysis 

To further investigate the mode lifetime for ^Hya, we com- 
pared the observed cumulative distribution of the power spec- 
trum from 1-190 /iHz with that from simulations (Fig.|9j. Note 
that the simulator is normalized so that on average it reproduces 
the observed power regardless of the mode lifetime (Paper 
I, Eq. 7). For each level in power, we calculated the frac- 
tion of t he power spectrum that is above t hat level (similar to 
Fig. 2 into elache & Scherrer 1983': see also lBrown et al.ll99l[ 
bedding et al-2005b: Bruntt et al..2005l) . We see only few time 
series of t = 17 days that have a power distribution similar to 
that observed and on average the difference is significant, while 
the distributions for r = 2 days resemble the observations much 
better. We note that there is a large intrinsic variation seen in the 
power distributions for r = 17 days. 

We also investigated the reliability of the observed frequen- 
cies (Table[3, using the simulations described in Sect. 13.21 The 
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Fig. 8. Distribution of the minimized scatter, crx, based on 10 
frequencies from each of 100 simulations with different mode 
lifetimes: 17 days, 5 days, 2 days, and 6 hours. The dashed 
line indicates the value from the observations, and the vertical 
arrows show the median. 



measured frequencies, together with the input frequencies and 
noise level are plotted in Fig. ^| Apart from the broadening 
of the mode frequencies due to damping, we also see false de- 
tections of alias peaks. This is most easily seen in panel (a) 
where the damping is less (a few examples of false detections 
are indicated). Our test shows that frequencies are not unam- 
biguous if S/N<7-8, even for mode lifetimes of 17 days. For 
short mode lifetimes it looks rather hopeless to measure the in- 
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Fig. 9. Panel (a): Cumulative distribution of power spectra 
from 100 simulations with r = 17 days (grey lines). The av- 
erage distribution is shown with a black line. The dashed line 
indicates the distribution from the observations. Panel (b): 
Similar to panel (a), but with r = 2 days. The noise level in 
the observed power spectrum is indicated. 
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dividual mode frequencies with useful accuracy. The observed 
frequencies are all with S/N<6.3 and hence cannot be claimed 
to be unambiguous. Hence the frequencies in Tablen(see also 
Fig.^Jbottom panel) should not be used for a direct compari- 
son with individual frequencies from a pulsation model. The 
lack of a clear comb pattern similar to Fig.^|(top panel) in the 
observed frequencies also supports a short mode lifetime. 



Fig. 10. Each panel shows 1000 measured frequencies and 
their amplitudes (10 from each of 100 simulations). Dashed 
lines indicate 7 times the input noise, and the solid grey lines 
are the input frequencies. Panel (a): Mode lifetime - 17 days. 
Panel (b): Mode Ufetime = 2 days. Panel (c): Mode lifetime - 
6 hours. Panel (d): Observed frequencies and the best match- 
ing comb pattern (grey lines). 



5. Discussion and future prospects 

Our best match for the mode lifetime of ^ Hya (t ~ 2 days) is 
in good agre ement with our estimate in Paper I but disagrees 
with theory iHoudek & Gou gh 2002). This discrepancy indi- 
cates that red giants could be used to better understand the 
mechanisms of the driving and damping of oscillations in a 
convective environment. The short mode lifetime of ^ Hya also 
suggests a narrow range of mode lifetimes for a large range 
of stars, from the main sequence to the red giants, and hence 
a steep decline in the 'quality' factor (Fig. The transition 
from the short mode lifetime regime to the much longer life- 
times of the semi-regular variables still needs to be investi- 



gated. It is likely to involve some kind of interaction between 
pure stochastic excitation and excitation by the k mechanism, 
which is respon sible for the oscillations we see in Mira stars 
(lBeddingl2003l) . 

A difficulty in using the current observations of ^Hya for 
asteroseismology arises from the severe crowding in the power 
spectrum. We now discuss possible origins of the crowding and 
some aspects of this issue. 

We do not expect crowding from high amplitude non-radial 
modes in the power s pectrum although low amplitude modes 
cannot be excluded dChristensen-Dalsgaard .20 04, Paper I). 
Additional simulations are needed to quantify the effect on the 
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Fig. 11. Power spectrum of simulated time series of ^ Hya 
(t - 2 days) using a window function from a two-site cam- 
paign where the sites are separated by 12 hours in longitude 
(see text). Dotted Unes indicate the input frequencies. 



lifetime estimate from low amplitude non-radial modes, but we 
expect it to be small. 

The crowding in the power spectrum comes partly from the 
single-site spectral window, which emphasizes the importance 
of using more continuous data sets from multi-site campaigns 
or space missions. To illustrate this, we construct the combined 
time series from a hypothetical two-site observing campaign on 
^Hya where the observing window from each site is identical 
to that of the present data set obtained with the Coralie spectro- 
graph. The other site is assumed to be a twin at the complemen- 
tary longitude (12 hour time shift), and the two individual data 
sets have identical sampling and noise. We simulated the stellar 
signal using t = 2 days, with other parameters the same as for 
Fig. 12 in Paper I. A power spectrum of such a two-site time se- 
ries is shown in Fig. ^2 Each mode profile is seen much more 
clearly, though slightly blended due to the short mode lifetime. 
Obviously, more can be obtained from such a spectrum than 
from our present data set (Fig. but a thorough analysis of 
similar simulations should be done to determine the prospects 
for doing asteroseismology on ^ Hya. 

As more high-quality asteroseismic data become available, 
simulations will continue to be an important tool for interpret- 
ing the data. However, one has to be careful about what can be 
deduced from simulations. Even though we may think, from an 
ideal perspective, that some parameters do not affect our mea- 
surements, this may turn out not to be true. A realistic noise 
source, the right window function, number- and scatter of fre- 
quencies, their amplitude and mode lifetime are all parameters 
that have to be considered when simulating oscillation data. 
Our comparison of the results from simplified simulations with 
more comprehensive simulations shows clearly that using sim- 
plified simulations can lead to erroneous conclusions. As an ex- 
ample, we measured an increase by a factor of two of min(crx) 
between simulations including only 1 and 18 input frequen- 
cies. This demonstrates that it is important to use realistic (non- 
simplified) simulations to relate frequency scatter to the mode 
lifetime. 



A complete application of our method on other stars (e.g. 
a Cen A, v Ind, t] Boo, and the Sun) will be published in a forth- 
coming paper. 

6. Conclusions 

We introduced a new method that uses only the scatter of 
the measured frequencies from a comb pattern to measure the 
mode lifetime of the red giant ^ Hya. This method takes into 
account false detections of noise and alias peaks. We find that 
the most likely mode lifetime is about 2 days, and we show that 
the theoretical prediction of 17 days dHoudek & Goughl2002h 
is unlikely to be the true value. 

Due to the high level of crowding in the power spec- 
trum, the signature from the p-modes is too weak to deter- 
mine the large separation to very high accuracy. However, 
our measurement sup ports the separation of 6.8 yuHz found by 
iFrandsen et"aD ( l2002l) . 

We conclude that the only quantities we can reliably ob- 
tain from the power spectrum of ^ Hya are the mode amplitude, 
mean mode lifetime, and the average large frequency separa- 
tion. 

Our simulations show that no ne of the m easured frequen- 
cies from the ^Hya data set (Frandsen et aL 200Z) can be 
regarded as unambiguous. Hence we suggest the measured 
frequencies given in Table are not used for direct mode- 
matching with a pulsation model to constrain the stellar model. 
Only in the case of a greatly improved window function could 
this be possible. 

Acknowledgements. This work was supported in part by the 
Australian Research Council. 

Appendix A: CLEAN test 

The method described in Sect. |3] required extraction of many 
frequencies. We therefore tested two different sine-wave fitting 
or CLEANing methods, simple CLEAN and CLEAN by simul- 
taneous fitting. We call them CLEANl and CLEAN2, respec- 
tively. The tests were done on our simulated time series to gain 
better understanding of the results from CLEANing and to de- 
termine which method was most favorable in our case. Both 
methods subtract one frequency at a time (the one with high- 
est amplitude), but CLEAN2 recalculates the parameters (am- 
plitude, phases, and frequencies) of the previously subtracted 
peaks while fixing the frequency of the latest extracted peak. In 
this way, the fit of the sinusoids to the time series is done si- 
multaneously for all peaks. CLEANl, however, does not recal- 
culate the parameters of previously subtracted peaks. The time 
used by CLEAN2 is a factor of 1.5(A^extract +1) longer than for 
CLEAN 1 , where A^extract is the number of extracted frequencies. 

We made simulations similar to those described in Sect. 13. 21 
but with 17 equally spaced frequencies and no noise added. 
From each set of 100 time series with a given mode lifetime, 10 
frequencies were extracted providing 1000 frequencies which 
we compared with the input frequencies. For coherent oscilla- 
tions CLEAN2 is doing a perfect job while CLEANl detected 
1% false peaks. If only true detections are considered, the scat- 
ter of the extracted frequencies relative to the input is roughly 
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equal to the frequency resolution (l/Tobs) for CLEAN2, while 
the CLEAN 1 results scatter twice as much. 

For non-coherent oscillations, the number of false detec- 
tions increases for both methods and the relative success rates 
of the two methods become equally good (or bad), for mode 
lifetimes in the order of the observational window or shorter 
In this regime, the scatter of the extracted frequencies gets 
dominated by the mode damping and is independent of the 
CLEANing method. For mode lifetimes below 5 days, the scat- 
ter of the frequencies is so pronounced relative to the spacing 
between modes and the alias peaks that they overlap and it can- 
not be determined whether an extracted frequency is false or 
not (see e.g. Fig.lsj"). 

The amplitude of a detected frequency found by both 
CLEANing methods generally differ by maximum +/-3Q% but 
can for a few cases differ by 70%, being most severe for shorter 
mode lifetimes. There is a general trend that CLEAN2 ascribes 
lower amplitudes to the peaks of highest amplitude in a spec- 
trum, and higher amplitudes to low amplitude peaks relative 
to CLEAN!. This is a result of the recalculation of the oscilla- 
tion parameters in CLEAN2. The detected amplitude difference 
was not considered as crucial as we were interested in using 
the frequencies only. Since there is observational evidence that 
the mode lifetime of ^Hya is below 5 days (Paper I), where 
the two methods perform equally well, we choose CLEAN 1 
due to its faster algorithm (by a factor of 17 in our case where 

A^extmct = 10). 

Appendix B: Robustness tests 

In this appendix we discuss each of the points (1-6) listed in 
Sect. l3.3l in turn: 

(1) In the minimization process, the initial guesses of Avq 
and Xq could give quite different outputs for Avq on a data set 
with a shallow dip in crx (i.e weak comb pattern). For data with 
a deep dip in crx (i.e. strong comb pattern), Avq is unaffected. 
We found that min(crx) was well determined if the initial guess 
of Avo was within ~ 1 //Hz of the true value. As a representative 
example we see ±1.5 //Hz change in Avq and only ±0.008 //Hz 
change in min{crx) for the set of observed frequencies, which 
produced a shallow dip in min(cr;f). Our final setup (plotted in 
Fig.|8|l was to fix the initial guess of Avq at 6.8 //Hz, which we 
believe is within 1 //Hz of the true value, and to use 7 different 
Xq values (separated 1 /(Hz apart). We then chose the solution 
with the lowest minimum found from these 7 trials. 

(2) We investigated the changes in the output parameters 
(min(crx), Avo), from varying the number of measured fre- 
quencies from 5 to 10. The average Avo is unaffected by the 
number of frequencies, but more frequencies give a more ac- 
curate Avo determination. For the observed frequencies Avq 
changes by ±0.5 //Hz in the tested regime (from 5 to 10 fre- 
quencies). Furthermore we observed higher values of min(crx), 
thus weaker comb patterns, with increasing number of frequen- 
cies, but the relative difference between the min(cr_y) distribu- 
tions from simulations and the observational data point did not 
change significantly. 

We decided to use 10 frequencies since we trusted our 
min(crx) distributions more when the Avo determinations were 



closer to the true value. Measuring more than 10 frequencies 
showed no clear advantage, for this particular data set. 

(3) Assigning weights to each frequency according to its 
S/N makes the absolute value of min(crx) slightly less sen- 
sitive to the number of measured frequencies used to calcu- 
late min(<Tx), which is expected. Due to the higher weight 
given to the central part of the damping profile for each mode 
dAnderson et al.ll99flh . the frequency separation is also slightly 
better determined using weights. However, mm(crx) will scatter 
slightly more, making it more difficult to distinguish min(crx) 
distributions from different mode lifetimes. This is probably 
a result of assigning lower weight to the tails of the damping 
profile. We chose not to use weights, to obtain the best possible 
result on the mode lifetime. 

(4) Changing the large separation of the input frequencies 
did not produce any change in mm(o-x) for the tested range 
(6.8-7.2 //Hz), but only changed the output Avq accordingly. 
Hence, it is not crucial for our results on the mode lifetime to 
know the true frequency separation to a very high accuracy. 

(5) We also changed the deviation, X(v,), from a comb pat- 
tern of the input frequencies. We used input frequencies that 
had a deviation twice as large as in the pulsation model (see 
Fig.|4}, and also a frequency set that followed a perfectly reg- 
ular comb. No significant change of the mm{o-x) distributions 
(Fig. |S} was observed. We believe this is because min(crx) is 
dominated by damping for the mode lifetimes we investigated. 

(6) Finally, we investigated whether the choice of frequency 
extraction method affects min(crx) and Avq. Using three dif- 
ferent algorithms on the observations provided 3 different sets 
of measured frequencies, and hence 3 different min(crx) val- 
ues that ranged from 0.23-0.28. The corresponding output Avq 
ranged from 6.7-7.2 //Hz. We note, that the min(crx) distribu- 
tions seen in Fig.|8]did not change significantly. 
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